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POLYNOMIAL APPROXIMATIONS OF THERMODYNAMIC PROPERTIES 
OF ARBITRARY GAS MIXTURES OVER WIDE PRESSURE 
AND DENSITY RANGES 

By Dennis O. Allison 
Langley Research Center 

SUMMARY 

Computer programs for flow fields around planetary entry vehicles require real-gas 
equilibrium thermodynamic properties in a simple form which can be evaluated quickly. 

To fill this need, polynomial approximations were found for thermodynamic properties of 
air and model planetary atmospheres. A coefficient-averaging technique was used for 
curve fitting in lieu of the usual least-squares method. The polynomials consist of terms 
up to the ninth degree in each of two variables (essentially pressure and density) including 
all cross terms. Four of these polynomials can be joined to cover, for example, a range 
of about 1000 to 11 000 K and 10 - 5 to 10® atmosphere (1 atm = 1.0133 x 10^ N/m2) for a 
given thermodynamic property. Relative errors of less than 1 percent are found over 
most of the applicable range. 


INTRODUCTION 

Computer programs used to analyze the flow field around a vehicle traveling at plan- 
etary entry velocities require real-gas equilibrium thermodynamic properties in a simple 
form which can be evaluated quickly. In the preliminary design studies for such vehicles, 
thermodynamic properties are needed for chemically reacting gas mixtures similar in 
complexity to air. However, the elemental compositions of other planetary atmospheres 
are not yet known, and the proposed models change as new discoveries are made. Further- 
more, the flow fields around blunt body shapes, which are being considered for entry vehi- 
cles, require the thermodynamic properties over wide ranges of pressure and density. A 
simple and versatile means of handling the thermodynamic properties consistent with the 
above considerations is needed. 

One approach is to read large tables of properties into the computer memory for 
interpolation purposes. However, the storage required by such tables can be prohibitive 
in view of the number of other variables which must be carried in the high-speed memory 
for a detailed flow-field calculation, especially in three-dimensional problems. Another 
approach is to represent the tables by two-dimensional approximate functions. (For such 



a problem, see ref. 1.) For air, Grabau (ref. 2) determined approximate functions by 
constructing two-dimensional continuous equations from straight-line segments, transi- 
tion functions, and eighth-degree polynomials. However, Grabau's approach requires 
tedious hand calculations which would have to be carried out for each model planetary 
atmosphere. Another approach, the one -dimensional spline fit (ref. 3), was used in ref- 
erence 4 but is not easily adaptable to two independent variables such as pressure and 
density (ref. 5). Also, in references 2 and 4 approximate functions are determined for 
first-order properties only and are differentiated to obtain second-order properties. 
(Second-order properties are those which depend on the partial derivatives of the species 
concentrations, for example, heat capacity and speed of sound. (See ref. 6.)) 

The present two-dimensional curve-fit procedure is designed to convert thermody- 
namic properties to a simple form suitable for automatic computation for various model 
planetary atmospheres. This procedure employs a 100-term polynomial of ninth degree 
in each of two variables to approximate any thermodynamic property. Since the input data 
(refs. 6 and 7) are not experimental, the type of smoothing of the least-squares method 
(which is discussed in ref. 8, for example) is not appropriate in determining the polynomial 
coefficients. That is, many thermodynamic properties are characterized by large varia- 
tions due primarily to chemical reactions; a least-squares curve fit tends to smooth out 
variations in computed data just as it smooths out scatter in experimental data. In an 
attempt to preserve the variations but to smooth out oscillations which may occur between 
input data points, the coefficients are determined by taking the average of polynomials 
which exactly reproduce alternate points of an input data array. No error criterion is met 
uniformly over the entire range of input data. Instead, large errors are concentrated on 
the edges of the range, and good accuracy is found throughout the central area. The cen- 
tral areas of several polynomials can be joined to accurately cover wide pressure and 
density ranges. Curve fits of the present type have been applied in the calculation of a 
blunt-body flow field in reference 9. 


SYMBOLS 

Because of the lack of lower case letters and subscripts on the computer, some of 
the symbols in the figures are different from those in the text. These symbols are indi- 
cated in parentheses in the symbol list. 


a (A) 

speed of sound 

a Q (AO) 

reference speed of sound, [1.4 — ] 

V V 

A n (x) 

polynomial coefficients (see eq. (2)) 


1/2 


2 



A n( x ) J A n( x ) 


coefficients whose average is A n (x) 


B m,n polynomial coefficients (see eq. (3)) 

Cp (CP) heat capacity at constant pressure per mole of undissociated mixture 
h (H) enthalpy of mixture per mole of undissociated mixture 

i (IX) index 

j (IW) index 


Mu 

P 

Po 

R 

s (S) 
T 


w 


molecular weight of undissociated mixture 
pressure 

reference pressure, 1.0133 x 10^ N/m2 (1 atm) 
universal gas constant, 8314.3 J/kmole-K 
entropy of mixture per mole of undissociated mixture 
absolute temperature 

reference temperature, 273.15 K 

W - Wi 1 
normalized W, 

W 2 o - Wi 2 


W = log 10 | 



Wj jth value of evenly spaced W's (see fig. 1) 


x 


normalized X, 


X - x x 
X 20 “ X 1 


1 

2 


X = log 



3 



ith value of evenly spaced X's (see fig. 1) 


Y = log 10 T 

Z any thermodynamic property 


J i,j 


input thermodynamic -property value corresponding to and Wj (see 
fig. 1) 


P 

Pr 


mass density 


reference mass density, 


M IjPo 


RT r 


CURVE -FIT PROCEDURE 

The analytical form of the curve -fit polynomial will be given, and will be followed 
by a discussion of some numerical studies used to determine parameters which give use- 
ful results. The problem is to find a procedure which gives accurate results over fairly 
wide ranges of two independent variables and which can be programed so that as model 
atmospheres are updated, the new versions can be treated automatically. 

The input thermodynamic properties can be obtained from a number of sources, and 
therefore, the calculation of these properties is not described herein. The properties for 
all the curve fits in this report were computed by use of the methods in references 6 and 7 

Analytical Form 

The basic curve-fit polynomial for any thermodynamic property Z is of ninth 
degree in each of two variables x and w, that is, 

9 9 

Z(x,w) = J ^ B m>n x%n ( 1 ) 

n=0 m=0 

The coefficients B m n are determined in two steps. First, for a fixed value of x, 

9 

Z(x,w) = ^ A n (x)w n (2) 

n=0 


4 



is solved for the 10 coefficients A n (x). Next, 


9 

A n (x)= £ B mn x m (3) 

m=0 

for n = 0, 1, . . 9 is solved for the 100 coefficients B mjn . In both of these steps, the 

coefficients are determined by taking the average of polynomials which exactly reproduce 
alternate points of an input data array. That is, the thermodynamic property Z(xi,Wj) 
is known for 20 evenly spaced values of Xj and wj. For a fixed value of x = xj, 

10 simultaneous equations are formed by substituting Z(xj,wj) and wj for 
j = 1, 3, . . .,19 into equation (2) and are solved for a set of 10 coefficients An(xj). 
Another set, An(xj), is determined for j = 2, 4, . . 20. A final set of A n (xj) is found 

by simply averaging the first two sets, 

AnW - ^4 - — (4) 


for n = 0, 1, . . .,9. The B m n are found in the same way from equation (3). The 
details of using the 400 input property values Zy in figure 1 to determine the 100 coef- 
ficients B mj n are given in the appendix. The relationship between the variables in 
equation (1) and those in figure 1 is 


X - Xj 

*20 - X T 




1 

2 


and 




W - Wi 

w = 

W 20 - W X 


1 

2 


( 5 ) 


In order to assess the accuracy of the fit whose form is given by equations (1) to (5), 
both relative and absolute errors are monitored for each property as follows: Relative 
errors (percent deviation from input values) are displayed by a chart of error codes, 
defined in table I. A chart (fig. 2, for example) contains one error-code entry for each 
input Zj j in figure 1. Absolute errors are checked by a series of six plots (fig. 3, for 
example) which compare curves from equation (1) with the input Zj j values. These 
curves cover the range of independent variables as follows: Two of the plots cross through 
the center and four pass near the edges of the X and W range. 

The large error codes found on the edges of a chart such as figure 2 are to be dis- 
regarded since the valid range is taken to be X 2 = X = Xjg and W 2 = W = Wjg, which 
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excludes the outermost rows and columns. Similarly, for a series of six plots such as 
figure 3, the deviation of each curve from the first and 20th input data points is to be dis- 
regarded. The valid range of each curve is taken to be 2 ~ i § 19 or 2 = j = 19. 

Greater accuracy could be achieved by choosing a smaller valid range. When the valid 
ranges of polynomials are joined as in figure 4, the 19th column of one chart lies on top of 
the second column of- the next chart. The maximum of the two error codes is given for 
each entry in that column. 


Numerical Studies 


Numerical studies show that the accuracy of the fit varies with choice of independent 
variables, degree of polynomial, coefficient averaging, and range covered by the indepen- 
dent variables. The present numerical studies were not exhaustive; therefore, the best 
combination of all these aspects may not have been found. The accuracy resulting from a 
good combination is illustrated in figures 4 and 5, which are for h/RT of air. The inde- 

'p\ /P/Po\ 


pendent variables are X = log^of^ - ) anc * W = logjof 


, the polynomial is of ninth 


;o/ \ p / p 0 ) 

degree, each coefficient is an average of two coefficients, and the ranges of the variables 
for a given polynomial are AX ~ 5.0 and AW ~ 0.4. Note that these four aspects are 
not entirely unrelated. For example, if the range of a variable is restricted, an accurate 
fit can be made with a lower order polynomial. However, in several subsequent examples, 
three of these aspects will be kept the same as given here while one will be changed. In 
each of these examples, the resulting error-code chart will be compared with that in fig- 
ure 4. 


Independent variables .- One alternate set of independent variables which was tried 
is X and Y = logjg T. Error codes for h/RT of air are given in figure 6. The X 
range is the same as that in figure 4, whereas the Y range is similar to that for W. 

A comparison of figure 6 with figure 4 shows that the errors in figure 6 are larger. This 
difference indicates that the independent variables X and Y are not as effective as X 
and W. 


Degree of polynomial .- The choice of the degree of the polynomials is arbitrary; 
ninth-degree polynomials give accuracy comparable to that of the input thermodynamic 
properties (refs. 6 and 7). The use of lower order polynomials might be desirable to save 
computer storage space, provided that sufficient accuracy is achieved for the intended 
application. Figure 7 gives error codes for seventh-degree polynomials approximating 
h/RT of air in the same X and W ranges as in figure 4. The maximum error code 
is 2 in both figures 7 and 4. However, more error codes of 1 and 2 appear for seventh- 
degree polynomials than for ninth-degree polynomials. 
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Coefficient averaging.- The coefficient averaging is eliminated by using a 10 by 10 
array of input data which consists of entries for every other X and every other W 
from the 20 by 20 array in figure 1. A ninth-degree polynomial exactly reproduces these 
100 input data points; these points must therefore have zeros for error codes in figure 8, 
which is for h/RT of air. The remaining 300 error codes are a measure of the oscil- 
lations between the 100 input data points. Figure 8 corresponds to the second polynomial 
of figure 4. It can be seen that the errors are larger (particularly for W = 0.9630) for 
the case without averaging. 

Range covered by a polynomial . - An attempt was made to approximate a thermody- 
namic property over very wide X and W ranges with one two-dimensional polynomial. 
An X range of Xjg - X 2 ~ 8 is covered in figure 9, compared with the previous range 
of Xjg - X 2 ~ 5 in figure 4. The range of W for the error codes of the polynomial in 
figure 9 is about the same as for the four polynomials combined in figure 4. The errors 
shown in figure 9 are judged to be much too large to justify the use of one polynomial to 
approximate a property over such wide ranges of variables. Although this case is 
extreme, it illustrates that the accuracy deteriorates rapidly as the range is increased. 

For applications in which higher pressures are encountered, variations in thermo- 
dynamic properties are generally weaker and are consequently easier to approximate. 
Smaller errors are shown in the valid range of figure 10, which is for -3 < X < 2, than in 
the third polynomial of figure 4, which is for -5 < X < 0. The W range is the same in 
both figures. 

DISCUSSION OF RESULTS 


The discussion will be limited to a few sample results for several thermodynamic 
properties of air and carbon dioxide. For all these results, the independent variables are 


x = lo sio(p] 

v 0/ 


and 


W = log 10 



the polynomial is of ninth degree, each coefficient 


is an average of two coefficients, and the ranges of the variables for a given polynomial 
are AX ~ 5.0 and AW ~ 0.4. Second-order properties are more difficult to approximate 
than first-order properties; however, some results will be given for both types. Polyno- 
mials approximating h/RT and c^/r are generally found to contain the largest errors 
of the approximations for all first- and second-order properties, respectively. Therefore, 
results for these two properties will be emphasized. 


Properties of Air 

Temperature .- The polynomial approximations of thermodynamic properties are 
functions of pressure and density but not temperature. If, for example, properties are 
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needed at a given temperature and density, the polynomial for temperature can be used in 
an iterative process to determine the corresponding pressure. All other properties can 
then be evaluated from the computed pressure and known density. This process is accu- 
rate because temperature is accurately approximated, as illustrated by the typical results 
in figures 2 and 3. (Recall that the outermost rows and columns of an error-code chart 
such as figure 2 are outside the valid range of the polynomial.) 

Entropy.- Entropy is also an easy property to fit. Figures 11 and 12 indicate that 
entropy is accurately approximated. 

Enthalpy .- The approximation of h/RT of air is shown in figure 4, where the valid 
ranges of four polynomials have been joined to cover a wide W range (corresponding to 
temperatures of about 1000 to 11 000 K). The edges of each error-code chart have been 
omitted from figure 4. Where the 19th column of one chart coincides with the second col- 
umn of the next chart, the maximum of each pair of error codes is shown. Note that most 
of the relative errors are less than 1 percent and none are larger than 2 percent in the 
valid range for h/RT. The plots in figure 5 correspond to the error-code chart of the 
third polynomial in figure 4, where the largest error codes are found. 

Heat capacity.- Second-order properties are more difficult to approximate because 
they have stronger variations. Results for one such property, c p y/R, are given in fig- 
ures 13 and 14. Figure 13 is constructed from error codes for the valid ranges of four 
polynomials in such a manner that the X and W ranges are the same as in figure 4. 

The plots corresponding to the third polynomial are given in figure 14 because the largest 
relative errors occur there. Note, for example, that the first row (X 2 = -5.1112) in fig- 
ure 13 contains a large error code of 5 where the value of c p ^R is small at IW = 12 in 
figure 14(a). In this case, the large error code is misleading since the absolute error is 
small. A similar situation exists for X^g = -2.7408 for IW = 17, 18, and 19. Most of 
the valid range contains relative errors of 1 percent or less; reasonably small absolute 
errors are found throughout the valid range. 

Speed of sound.- The X and W range for the a/ao results in figures 15 and 16 
is the same as that of the third polynomial in figure 13 because the largest errors in a/ao 
are found there. Most of these errors are less than 1 percent, and none are larger than 
5 percent in the valid range. 


Properties of CO 2 

Results similar to those for air were found for the model Mars atmospheres con- 
sidered in reference 9, where the present method was used. Results for those mixtures 
of carbon dioxide (CO 2 ), oxygen (O 2 ), and argon can be illustrated by considering 100 per- 
cent CO 2 . Actually, the dissociation of 100 percent CO 2 is more complicated (hence, 
curve fitting is more difficult) than that of air. This complexity is caused by the three 


8 



major dissociation processes for CO 2 (2CC>2 — 2CO + O 2 , O 2 — 20, and CO — C + O) 
compared with the two for air (O 2 — 20 and N 2 — 2N). 

Enthalpy . - Figures 17 and 18 are for h/RT of CO 2 , which is negative at tempera- 
tures below about 3000 K. The negative values pose no problem in curve fitting except 
that the error code is meaningless where h/RT passes through zero. This problem has 
been avoided in the error codes of figure 17 by simply adding a constant to h/RT to pre- 
vent the property from passing through zero. The first polynomial approximates 
h/RT + 45; the second, h/RT + 15; and the last two, h/RT. The X and W ranges in 
figure 17 are the same as those for air in figure 4. The largest of the relative errors in 
the valid range is 4 percent, and most of the errors are less than 2 percent. Since the 
largest error codes occur in the third polynomial, plots for it are given in figure 18. 

Heat capacity .- Results for Cp^R of CO 2 are given in figures 19 and 20 for the 
same X and W ranges as before. In figure 19, as in previous figures, the largest 
error codes occur in the third polynomial, for which plots are given in figure 20. These 
plots indicate that the value of Cp/R is small in areas where the error codes are largest. 
In these areas the plots in figure 20 show that the absolute errors are reasonably small. 

CONCLUDING REMARKS 

Ninth-degree two-dimensional polynomials can be used with a coefficient-averaging 
technique to approximate thermodynamic properties of air and model planetary atmo- 
spheres. For a temperature range of about 1000 to 11 000 K and a pressure range of 10“5 
to 10° atmosphere (1 atm = 1.0133 x 10^ N/m2), four polynomials can be joined to approx- 
imate a single thermodynamic property. These four polynomials would be stored in a 
computer as 400 polynomial coefficients. Relative errors of less than 1 percent are found 
over most of the applicable range. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 27, 1972. 
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APPENDIX 


DETERMINATION OF POLYNOMIAL COEFFICIENTS 

The 400 values of any property Z, as illustrated in figure 1, are used to determine 
the 100 coefficients of the two-dimensional ninth-degree polynomial 


Z(x,w) 



B 0,0 

B 0,l 

b 0,2 

• • • b 0,9 

1 


B l,0 

b m 

b 1,2 

■ ; • b 1,9 

w 

99 r 1 

^ ^ B m>n x m w n =1 x x2 ... x9 

m=0 n=0 

b 2,0 

b 2,1 

b 2,2 

• • • b 2,9 

w2 


b 9,0 

b 9,1 

b 9,2 

• • • b 9,9 

w 9 


where 


and 


X - Xj 
X 20‘ X 1 


1 

2 




W - Wj[ j 
W = W 2 o- w i-2 


(Al) 


(A2) 


These variables x and w are normalized to range between -0.5 and 0.5 to avoid the 
possibility of ill- conditioning (ref. 10) in matrices whose elements are powers of x's or 
w's. The normalization is also numerically convenient, since the x's and w's always 
have the same values: xj = Wj = -0.5, X 2 = W 2 = -0.44737, Xg = Wg = -0.39474, . . ., 
x 18 = w 18 = 0.39474, xjg = w^g = 0.44737, and Xgg = Wgg = 0.5. 

The solution for the 100 coefficients B m n of equation (Al) begins by determining 
two "exact" one-dimensional polynomials for each given x value, with w as the vari- 
able. (An "exact" polynomial is referred to as a polynomial interpolating function in 
chapter 3 of reference 11. It passes through every point.) One polynomial is required to 
exactly reproduce the first, third, . . ., and the 19th property values while the other exactly 
reproduces the second, fourth, . . ., and 20th values. For the ith x value, one of these 
polynomials is found by solving 
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APPENDIX - Continued 


Z M 

Z l,3 

Z l,5 


Z l,19 


a l,0 

a l,l 

a l,2 - * • 

a l,9 

1 

1 

1 

1 

z 2,l 

z 2,3 

Z 2,5 


Z 2,19 


a 2,0 

a 2,X 

a 2,2 

a 2,9 

W 1 

w 3 

w 5 . 

• w 19 

z 3,l 

Z 3,3 

Z 3,5 


Z 3,19 


a 3,0 

a 3,l 

a 3,2 • • • 

a 3,9 

WX 2 

W 3 2 

w g 2 . 

• w 19 2 

Z 4,l 

Z 4,3 

Z 4,5 


Z 4,19 


a 4,0 

a 4,l 

a 4,2 • • ’ 

a 4,9 





z 5,l 

Z 5,3 

Z 5,5 


Z 5,19 

_ 

a 5,0 

a 5,l 

a 5,2 - - - 

a 5,9 





z 6,l 

Z 6,3 

Z 6,5 


Z 6,19 


a 6,0 

a 6,l 

a 6,2 - • • 

a 6,9 


W 3 9 

w 5 9 . 

• w 19 9 

z 20,l 

Z 20,3 

Z 20,5 

. . . 

Z 20,19 


a 20,0 

a 20,l 

a 20,2 • • • 

a 20,9 





for the ith 

row of 

H,n 

and the other by solving 






z 1,2 

z l,4 

z l,6 


z l,20 


b l,0 

b l,l 

b l,2 • • • 

b l,9 

1 

1 

1 

1 

z 2,2 

z 2,4 

z 2,6 


z 2,20 


b 2,0 

b 2,l 

b 2,2 • • • 

b 2,9 

w 2 

w 4 

W 6 . 

• w 20 

Z 3,2 

z 3,4 

Z 3,6 


Z 3,20 


b 3,0 

b 3,l 

b 3,2 • • • 

b 3,9 

w 2 2 

w 4 2 

Wg2 . 

• w 20 2 

Z 4,2 

Z 4,4 

Z 4,6 


Z 4,20 


b 4,0 

b 4,l 

b 4,2 • • ’ 

b 4,9 

• 




Z 5,2 

Z 5,4 

Z 5,6 


Z 5,20 

= 

b 5,0 

b 5,l 

b 5,2 • • • 

b 5,9 





Z 6,2 

Z 6,4 

Z 6,6 


Z 6,20 


b 6,0 

b 6,l 

b 6,2 • • • 

b 6,9 

w 2 9 

w 4 9 

Wg 9 . 

9 

• w 20 

z 20,2 

z 20,4 

z 20,6 


z 20,20 


b 20,0 

b 20,l 

b 20,2 • • ■ 

b 20,9 






(A3) 


(A4) 


for the ith row of t>i >n . The left-hand side of equation (A3) is formed from the odd- 
numbered columns of figure 1; the left-hand side of equation (A4) is formed from the even- 
numbered columns. The two polynomials are then averaged to obtain one polynomial given 
by the ith row of equation (A5), 


Z(x 1; w) 


A l,0 

A M 

A l,2 

• • • A l,9 

1 

Z(x 2 ,w) 


A 2,0 

A 2,l 

A 2,2 

• * • A 2,9 

w 

Z(x 3 ,w) 


A 3,0 

A 3,l 

A 3,2 

* * • A 3,9 

w 2 

Z(x 4 ,w) 


A 4,0 

A 4,l 

A 4,2 

* « * A 4,9 

• 

Z(x 5 ,w) 

— 

a 5,0 

A 5,l 

A 5,2 

• * • A 5,9 

• 

Z(x 6 ,w) 


a 6,0 

A 6,l 

A 6,2 

• • • a 6,9 

w 9 

z (x20» w ) 


A 20,0 

A 20,l 

A 20,2 

• • • A 20,9 



where A i n = 31,11 — — for i = 1, 2, . . .,20 and n = 0, 1, . . .,9. 
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APPENDIX - Continued 


At this point, the ith row of j is approximated by a polynomial such that 
9 


Z( Xi ,w) - £ A i n w n 


n=0 

for i = 1, 2, 3, . . .,20. Now the nth column of Aj n is to be approximated by a polyno- 
mial function of x. Two "exact" polynomials are formed for the nth column of Aj n by 
first solving 



A l,0 

A l,l 

a 1,2 

. . . 

A 

1,9 


1 

x r 

x l 2 • 

• x l 9_ 

c 0,0 

c 0,l 

c 0,2 • 

• c 0,9 


A 3,0 

a 3,1 

a 3,2 


a 3,9 


1 

x 3 

x 3 2 

• x 3 9 

c l,0 

c l,l 

c l,2 • 

c l,9 


A 5,0 

a 5,1 

A 5,2 

* * * 

A 

5,9 

= 

1 

x 5 

x 5 2 

• x 5 9 

c 2,0 

c 2,l 

c 2,2 • 

■ c 2,9 


A 19,0 

A 19, 1 

A 19,2 


A 

19,9 


1 

x 19 

x 19 2 • 

x 9 
x 19 

c 9,0 

c 9,l 

c 9,2 • • 

c 9,9 

for the nth column of 

c m,n 

and then 

solving 






a 2,0 

a 2,1 

a 2,2 

. . . 

A 2,9 


1 

x 2 


V 

x 2 3 . 

x 2 9 

d 0,0 

d 0,l 

d 0,2 ' 

• d 0,9 

a 4,0 

A 4,l 

A 4,2 


A 4,9 


1 

X 4 

x 4 2 

x 4 3 . 

. . x 4 9 

d l,0 

d l,l 

d l,2 • 

• d l,9 

a 6,0 

A 6,l 

A 6,2 

* * * 

A 6,9 

= 

1 

x 6 

x 6 2 

x 6 3 . 

. . x 6 9 

d 2,0 

d 2,l 

d 2,2 • 

• d 2,9 

a 20,0 

A 20,l 

a 20,2 

... 

A 20,9 


1 

x 20 ' 

x 20 2 

x 20 3 • 

• • x 20 9 

d 9,0 

d 9,l 

d 9,2 • 

• d 9,9 


(A6) 


(A7) 


for the nth column of d m>n . The left-hand side of equation (A6) is formed from the odd- 
numbered rows of Aj n of equation (A5); the left-hand side of equation (A7) is formed 
from the even-numbered rows. The two polynomials are then averaged to obtain one 
polynomial given by the nth column of equation (A8), 


A 0 (x) A x (x) A 2 (x) 


A 9 (x) 


u 


1 x 


B 0,0 

B 0,l 

b 0,2 

• • b 0,9 

B l,0 

B l,l 

b 1,2 

• • b 1,9 

b 2,0 

b 2,1 

b 2,2 • 

• • b 2,9 

b 9,0 

b 9,1 

b 9,2 ■ • 

• b 9,9 


(A8) 
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APPENDIX — Continued 


where 


Bm,n 


c m,n + ^m,n 
2 


(A9) 


for m, n = 0, 1, . . ., 9. The two-dimensional polynomial approximation is finally 
formed by substituting B mjI1 of equation (A9) into equation (Al). 

The set of B mj n and Xj, X 20 , Wj, and W 20 for the third polynomial of fig- 
ure 4 for h/RT of air are given here, and a sample evaluation of the polynomial is given. 
The B m>n are ordered as follows: 


B 0,0 

B 0,l 

B 0,2 

B 0,3 

B 0,4 

B 0,5 

B 0,6 

b 0,7 

b 0,8 

b 0,9 

b 9,5 

b 9,6 

b 9,7 

b 9,8 

b 9,9 


The integer following the E is a power of 10; for example, 1.2081E+03 = 1.2081 x 10^. 


m = 0: 

1.8842E+01 

1.8724E+01 

-1.9274E+01 

-7.4597E+01 

-1.4161E+02 


2.2529E+01 

1.2081E+03 

1.8045E+03 

-2.4644E+03 

-5.0926E+03 

m = 1: 

-1.9223E+01 

5.9860E+00 

1.4767E+02 

4.5817E+02 

1.2088E+02 


-4.9694E+03 

-8.4302E+03 

1.5751E+04 

2.7259E+04 

-1.3093E+04 

m = 2: 

5.7650E-01 

-7.9885E+01 

-4.7190E+02 

-1.2060E+03 

7.5502E+03 


4.0263E+04 

-3.0717E+04 

-2.8412E+05 

2.8112E+04 

6.098 1E+05 

m = 3: 

1.5838E+01 

1.7215E+02 

7.5910E+02 

-3.2372E+03 

-3.7724E+04 


-3.8374E+04 

3.4992E+05 

6.2886E+05 

-9.1143E+05 

-1.8723E+06 

m = 4: 

-2.3310E+01 

-2.6170E+02 

1.8817E+03 

3.1578E+04 

1.0701E+04 


-5.8825E+05 

-4.4625E+05 

3.6933E+06 

1.7799E+06 

-7.4600E+06 

m = 5: 

1.8278E+01 

-2.6375E+02 

- 1.023 1E+04 

-4.6790E+04 

2.9272E+05 


1.5589E+06 

-2.3791E+06 

-1.3535E+07 

5.7175E+06 

3.4212E+07 

m = 6: 

3.1957E+01 

1.7667E+03 

5.066 1E+03 

-1.4085E+05 

-4.1131E+05 


2.3622E+06 

5.3044E+06 

-1.3483E+07 

-1.7221E+07 

2.4399E+07 

m = 7: 

-5.9066E+01 

-1.0437E+03 

2.9512E+04 

3.4704E+05 

-6.6843E+05 


-9.2474E+06 

4.2249E+06 

7.5263E+07 

-7.1884E+06 

-1.8508E+08 

m = 8: 

1.1476E+01 

-2.6431E+03 

-2.1784E+04 

1.9858E+05 

1.1113E+06 


-2.9386E+06 

-1.3065E+07 

1.3530E+07 

4.1034E+07 

-1.6564E+07 

m = 9: 

-2.2317E+01 

2.1902E+03 

-2.5838E+04 

-6.5313E+05 

2.5673E+05 


1.6827E+07 

1.6364E+06 

-1.3486E+08 

-1.3248E+07 

3.2885E+08 

For these values of B mjI1 , the Xj, X 20 ) 

Wi, and W 20 

are 



X 1 = -5.4075 

X 20 

= 0.2222 




Wi = 1.3310 

W 20 

= 1.7680 
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APPENDIX — Concluded 


This polynomial for h/RT can only be used in the valid range as shown in figure 4, 
-5.1112^X^-0.0741 and 1.3540 = W = 1.7450. If, for example, h/RT is desired at 
the point X = -3.0 and W = 1.5, equation (A2) is used to find that x = -0.07236 and 
w = -0.11327. These normalized variables and B mjll are then substituted into equa- 
tion (Al) to evaluate h/RT = Z(x,w) = 17.9. 
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Figure 1.- Array of input values of a 
thermodynamic property Z. 
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Figure 2.- Error-code chart (table I) for polynomial 
approximating T of air. 
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Figure 4.- Error codes for valid ranges of four polynomials joined to approximate h/RT of air. 


CD 



24 


O input data 
Polynomial 




I I I I i I I i i 

1 6 8 10 12 14 16 18 20 

IW 

(B) IX □ 10 (X □ -2.7408). 


L. I I I I I I | 

0 2 4 6 8 10 ' 12 14 16 

IX 

(E) IW = 10 (W □ 1 .5380) . 




IW 


(C) IX - 13 IX = -0.0741). (F) IW □ 13 (W □ 1.7450). 

Figure 5.- Comparison of third polynomial of figure 4 approximating 
h/RT of air with input data. 
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Figure 6.- Error codes for valid X and Y ranges of four polynomials joined to approximate h/RT of air. 
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Figure 7.- Error codes for valid ranges of four seventh-degree polynomials joined to approximate h/RT of air 
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Figure 8.- Error-code chart for polynomial without averaging 
approximating h/RT of air in same range as second 
polynomial of figure 4. 
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Figure 9.- Error-code chart for polynomial 
approximating h/RT of air over very 
wide X and W range. 
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Figure 10.- Error-code chart for polynomial approximating 
h/RT of air in same W range as that of third polyno- 
mial of figure 4 but in a higher X range. 
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Figure 11.- Error-code chart for polynomial 
approximating s/R of air. 
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Figure 12.- Comparison of polynomial approximating s/R of air with input data, 
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Figure 13.- Error codes for valid ranges of four polynomials joined to approximate Cp/R of air. 
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Figure 16.- Comparison of polynomial approximating a/a 0 
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Figure 17.- Error codes for valid ranges of four polynomials joined to approximate h/RT + Constant of CO2 
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Figure 19.- Error codes for valid ranges of four polynomials joined to approximate Cp/R of CC> 2 . 
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Figure 20.- Comparison of third polynomial of figure 19 
approximating c p /R of CO2 with input data. 
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